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We present meson-meson (Wilson loop) correlators in Z{2) center vortex models for the infrared 
sector of Yang-Mills theory, i.e., a hypercubic lattice model of random vortex surfaces and a con¬ 
tinuous 2+1 dimensional model of random vortex lines. In particular we calculate quadratic and 
circular Wilson loop correlators in the two models respectively and observe that their expectation 
values follow the area law and show string breaking behavior. Further we calculate the catenary 
solution for the two cases and try to find indications for minimal surface behavior or string surface 
tension leading to string constriction. 
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I. INTRODUCTION 

Many physical principles can be formulated in terms of variational problems. For example the least-action principle 
is an assertion about the nature of motion that provides an alternative approach to mechanics completely independent 
of Newton’s laws. Not only does the least-action principle offer a means of formulating classical mechanics that is 
more flexible and powerful than Newtonian mechanics, but also variations on the least-action principle have proved 
useful in general relativity theory, quantum held theory, and particle physics. As a result, this principle lies at the 
core of much of contemporary theoretical physics. 

Calculus of variations seeks to hnd the path, curve, surface, etc., for which a given function has a stationary value, 
which, in physical problems, is usually a minimum or maximum. In this article we want to use calculus of variations 
to analyze the (minimal) area law behavior of quark conhnement in meson-meson correlators. For this behavior the 
non-Abelian behavior of the gauge group is essential, and especially the center vortex model of conhnement has proved 
to be very successful in predicting the area law. 

Center vortices [H-Q are topological defects associated with the elementary center degree of freedom of the QCD 
gauge held. In D-dimensional space-time, they are (thickened) (D — 2)-dimensional chromo-inagnetic hux degrees of 
freedom and have a clear theoretical link to conhnement [I|, 0^3 • This was conhrmed by a multitude of numerical 
simulations, in lattice Yang:-Mills theory, see em. EMI , and within the infrared effective models of center vortices 
under investigation E3 t 23- Recent results have also suggested that the center vortex model of conhnement is 
more consistent with lattice results than other currently available models. Lattice QCD simulations indicate further 
that vortices are responsible for the spontaneous breaking of chiral symmetry (ySB), dynamical mass generation and 
the axial Ua{^) anomaly (23 - 1^ . and thus successfully explain the non-perturbative phenomena which characterize 
the infrared sector of strong interaction physics. 

In the present work we measure meson-meson correlators in effective random center vortex models, which will be 
introduced in the next section. Meson and baryon (Polyakov line) correlators were analyzed in the random vortex 
world-surface model in [d^, showing that the correlators follow an area law behavior. Here we investigate quadratic 
and circular Wilson (not Polyakov) loop correlators, focusing on the minimal area law. Motivated by the minimal 
surface of revolution problem we want to analyze possible signs of catenary solutions. 

The hanging chain or catenary problem (the world catenary comes from the Latin word catena meaning chain) was 
first posed in the Acta Eruditorium in May 1690 by Jacob Bernoulli as follows: To find the curve assumed by a loose 
string hung freely from two fixed points. Earlier, Galileo mistakenly conjectured that the curve was a parabola. Later 
Joachim Jung proved that the curve cannot be a parabola but without presenting any solution of the real curve. In 
June 1691 there were three solutions published, from Leibniz, Huygens and Johann Bernoulli brother of Jacob. Even 
though these mathematicians approached this problem in three different ways they concluded that the curve was the 
hyperbolic cosine, which then came to be known as the catenary. 

The catenoid is a surface of revolution with minimum area: Any surface of revolution generated by a different curve 
joining the same endpoints and having the same length, must have a larger surface area. The proof of this fact is less 
elementary and involves calculus of variations (see e.g. Hi]). In section uni we derive the so-called catenary solution 
using the calculus of variations, for both, standard circular and quadratic Wilson loops and compare the results with 
numerical measurements in section HVl We finish with concluding remarks inIVl 


II. RANDOM CENTER VORTEX MODELS 

Center vortices are expected to behave as random lines (for D = 3) or random surfaces (for D = A). The magnetic 
flux carried by the vortices is quantized in units which are singled out by the topology of the gauge group, such that 
the flux is stable against small local fluctuations of the gauge helds. The vortex model of confinement states that the 
deconfinement transition is simply a percolation transition of these chromo-magnetic flux degrees of freedom. The 
following infrared effective models in D = 3 and D = A dimensions have confirmed these expectations and in the 
latter case also account for topological susceptibility and the (quenched) chiral condensate. 


A. Random vortex world-line model 

In [ 21 I, [23 we introduced a model of random flux lines in D = 2 -|- 1 space-time dimensions. The lines are composed 
of straight segments connecting nodes randomly distributed in continuous three-dimensional space. Each node is 
connected to two lines, hence the conhguration only consists of closed vortex clusters. The physical space in which the 
vortex lines are defined is a cuboid x Lt with ’’spatial” extent Lg, ’’temporal” extent Lt and periodic boundary 
conditions in all directions. Within this paper we use volumes with Ls = Lt = 16, where hnite size effects are under 
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control. The vortex length L between two nodes is restricted to a certain range Lmm < L < Lmax- This range in some 
sense sets the scale of the model; for practical reasons we choose a scale of L Ri 1, i.e. Tmin = 0.3 and L^nax = 1-7 in 
these dimensionless units. All updates resulting in vortex lengths L out of the range Lmm < L < Lmax are rejected. 

An ensemble is generated by Monte Carlo methods, starting with a random initial configuration. Allowance is 
made for nodes moving as well as being added or deleted from the configurations during Monte Carlo updates. A 
Metropolis algorithm is applied to all updates using the action S = aL + with action parameters a = 0.11 and 
7 = 0.33 for the vortex length L and the vortex angle (p between two adjacent segments respectively. The difference of 
the action of the affected nodes before and after the update determines the probability of the update to be accepted. 
The move update moves the current node by a random vector of maximal length Vm = ^Lmin, it affects the action 
of three nodes, the node itself and its neighbors. The add update adds a node at a random position within a radius 
To = SLmin around the midpoint between the current and the next node. The action before the update is given by 
the sum of the action at the current and the next node, while the action after the update is the sum of the action 
at the current, the new and the next node. Deleting the current node, on the other hand, affects three nodes before 
the update and only two nodes after the update. Therefore the probability for the add update is in general much 
smaller than for the delete update and the vortex structure would soon vanish if both updates were tried equally 
often. Hence, the update strategy is randomized to move a node in two out of three cases (66%), and apply the add 
update about five times more often than the delete update (28% vs. 6%). A density parameter d is restricting the 
number of nodes in a certain volume. The add update is rejected if the number of nodes within a 3 x 3 x 3 volume 
around the new node exceeds the density parameter d. 

Furthermore, Monte Carlo updates disconnecting and fusing vortex lines were implemented, i.e., when two vortices 
approach each other, they can reconnect or separate at a bottleneck. The ensemble therefore will contain not a fixed, 
but a variable number of closed vortex lines or ’’vortex clusters”. Given that the deconfining phase transition is a 
percolation transition, such processes play a crucial role in the vortex picture. If the current node is not deleted, all 
nodes around the current node are considered for reconnections. The reconnection update causes the cancellation of 
two close, nearly parallel vortex lines and reconnection of the involved nodes with new vortex lines. The reconnection 
update is also subjected to the Metropolis algorithm, considering the action of the four nodes involved. 

The model was shown to exhibit both a low-temperature confining phase and a high-temperature deconfined phase, 
as well as phase transitions with respect to vortex density and segment length L. The predictions of the model for the 
spatial string tension in the deconfined phase quantitatively match corresponding SU (2) lattice Yang-Mills results. 


B. Random vortex world-surface model 

The random vortex world-surface model was introduced in [I^, |4^ based on the notion that the Yang-Mills vacuum 
is populated by collective magnetic vortex degrees of freedom which are represented by closed two-dimensional world- 
surfaces in four-dimensional (Euclidean) space-time. The chromomagnetic flux carried by the vortices is quantized 
according to the center of the gauge group. The vortex world-surfaces are treated as random surfaces, an ensemble 
of which in practice is generated using Monte Carlo methods on a hypercubic lattice; the surfaces are composed 
of elementary squares (plaquettes) on that lattice. The spacing of the lattice is a fixed physical quantity, related 
to an intrinsic thickness of the vortex fluxes, and represents the ultraviolet cutoff inherent in any infrared effective 
framework. The action governing the ensemble is related to the surface curvature: If two elementary squares which 
are part of a vortex surface share a lattice link but do not lie in the same plane, this costs an action increment c, where 
the parameter is set to c = 0.24, such as to reproduce the SU{2) Yang-Mills ratio of the deconfinement temperature 
to the square root of the zero-temperature string tension, Td^/^ = 0.69 [T^ . 

Physically, the random vortex surfaces represent quantized chromomagnetic flux. This means that they contribute 
in a characteristic way to Wilson loops; if one chooses an area spanning a given Wilson loop - the choice of area is 
immaterial due to the continuity of flux - then for each time a vortex world-surface pierces that area, the Wilson 
loop acquires a phase factor corresponding to the center of the gauge group. For the case of SU (2) color treated in 
this paper, the Wilson loop picks up a (—1) = exp(±i7r). Note that Wilson loops in the 4D model are defined on a 
lattice dual to the one on which vortices are defined, i.e. on a lattice shifted by the vector (a/2, a/2, a/2, a/2), where 
a denotes the lattice spacing. Thus, the notion of a vortex piercing a Wilson loop area is unambiguous. 

The vortex ensemble is generated subject to the constraint of continuity of flux (modulo 27r, ie., modulo Dirac 
strings forcing thee vortex surfaces to be closed, as already mentioned above. As described in detail 

in [20|, continuity of flux in practice is guaranteed during the generation of the vortex world-surface ensemble by 
performing updates simultaneously on the six squares making up the surface of an elementary three-dimensional cube 
in the lattice. This is done in a way which corresponds to superimposing the (continuous) flux of a vortex of the 
shape of the elementary cube surface onto the flux previously present. 

For the SU (2) gauge group, the model was shown to exhibit both a low-temperature confining phase as well as 
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a high-temperature deconfined phase [T^, separated by a second-order phase transition [2^. The predictions of the 
model for the spatial string tension in the deconfined phase , the topological susceptibility [i^ and the (quenched) 
chiral condensate quantitatively match corresponding SU(2) lattice Yang-Mills results. Building on this initial 
progress, the random vortex world-surface model was extended to the case of SU{3) color. Studies of the confinement 
properties yielded a weakly first-order deconfinement phase transition [^ . also seen in the vortex free energy (^ . 
which represents an alternative order parameter for confinement. Furthermore, a Y law for the baryonic static 
potential was observed (48l| and topological susceptibility was studied in (ssj j. which is instrumental in determining, 
via the C/y 4 (l) anomaly, the mass of the rj' meson. Spurred by investigations of Yang-Mills theories with a wider variety 
of gauge groups, aiming at a better understanding of possible confinement mechanisms (d^ . the confinement properties 
of the random vortex world-surface model were subsequently also studied for S't/(4) color and Sp{2) color M- 
These studies showed that the vortex picture can accomodate such diverse color symmetries, while indicating the 
limitations of the very simple effective dynamics which had proven adequate in the SU{2) and SU{3) cases. Finally, 
as already mentioned above, meson and baryon correlators were analyzed in the random vortex world-surface model 
for the latter cases in [d^, showing that the correlators follow an area law behavior and are in qualitative agreement 
with lattice studies. 


III. THE CATENARY SOLUTION 


A power line hanging between two poles shows us a curve called the catenary. The shape of this natural curve 
can be derived from a differential equation describing the physics behind a uniform flexible chain hanging by its own 
weight. The problem was also proposed by W. Symmond in [52| in 1894. J. C. Nagel gave a solution in that 
is easy to describe if we accept the reasonable fact that the catenary is the curve with the lowest center of gravity. 
The problem can be formulated in the same way as minimizing the surface of revolution of a curve between two fixed 
points, well known from calculus of variation, see e.g. (53 |: 

Given Pi, P 2 in the plane, find a curve y{x) from Pi to P 2 such that the surface of revolution obtained by revolving 
the curve about the cc-axis has minimum surface area. In other words, minimize 27r yds with y{xi) = Pi and ds 

the line segment along the curve y{x). If Pi and P 2 are not too far apart, relative to X 2 — xi then the solution is a 
Catenary (the resulting surface is called a Catenoid), see Fig.[T] 



FIG. 1. Solutions of the catenary problem for circular Wilson loops: a) Catenoid and b) Goldschmidt’s discontinuous solution. 


Otherwise the solution is Goldschmidt’s discontinuous solution (discovered in 1831) obtained by revolving the curve 
which is the union of three lines: the vertical line from Pi to the point (a:i, 0), the vertical line from P 2 to (x 2 , 0) and 
the segment of the x-axis from (xi,0) to (x2,0). This example also straightforwardly translates to the case of the 
minimal area law for the Wilson loops, i.e., if the two Wilson loops are far apart, we only measure the product of their 
expectation values or the string tensions of the two individual meson pairs. However, if the meson pairs approach 
each other, we measure mixed states all four mesons interact with each other and at some point the strings will flip 
and form new bound states. The question now is if we can interpret these mixed state measurements with catenary 
solutions. Therefore we derive the equations for the set of boundary conditions of our problem, starting with the 
circular Wilson loop. 
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A. Circular Wilson loops 


For circular Wilson loops the surface of revolution problem is straightforward. We want to find a curve from a point 
(—£)/2, R) to a point (I?/2, R), where R is the radius of the Wilson loop and D is the distance between two Wilson 
loops, which, when revolved around the x-axis yields a surface of smallest surface area A, i.e., the minimal surface. 
The area element is 


so the surface area is 


dA = 2'Kyds = 2'Ky\/l + y'‘^dx, 

rD/2 _ rD/2 _ 

A = 2tt y\/l + y^dx = An yy/l + y^dx, 

J-D /2 Jo 


and the quantity to minimize is 


/ = + 2/'^. 

This can be solved analytically by using the Beltrami identity, see appendix!^ and gives the solution 

y = acosh(^—^), 
a 

with two constants a and b to be determined from the boundary conditions 

„ ^-Dl2-b^ ,,Dl2-b^ 

R = acosh(-) = acosh(-). 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


Since cosh(—x) = cosh(x) we get D/2 — b = D/2 + b, hence & = 0, as it must by symmetry, and the minimal surface 
solution reduces to 


y = a cosh.(x/a) with the boundary condition R = a cosh{D/2a), 


( 6 ) 


to determine the constant a. But for certain values R and D, this equation has no solution. The mathematical 
interpretation of this fact is that the surface breaks and forms circular disks in each ring to minimize area, i.e. the 
Goldschmidt solution. Physically it simply means that the minimal area is given by the two Wilson loops rather than 
any surface of revolution, the two interpretations however do not agree exactly. In appendix |B] we derive that we 
only obtain catenary solutions for D/R < 1.325 (IB3I) . but not all solutions are absolute minima of our minimal area 
problem. The surface area of the minimal catenoid is given by © 


. 0/2 - 

A = An yyjl-\-y'^dji 
Jo 


but since according to Eqs. dH) and 


y = acosh(x/a) = \/l + y'^a, 


( 7 ) 

( 8 ) 


, 47r 

A= — 

a . 


1-0/2 pD/2 

/ y^dx = 47ra / coslJ^ {x/a)dx 

0 Jo 


. , ,2x. 2x 

smh( —) H- 

a a 


D/2 


= na^ [sinh(£)/a) -\- D/a] . 


J 0 


( 9 ) 

( 10 ) 


The surface area of the catenoid equals that of the Goldschmidt solution when (flUl) equals the area of two disks, the 
circular Wilson loops with radius i?, i.e.. 


naJ [sinh(Z)/a) + D/a] = 2nR^ (11) 

sinh(Zl/2a) cosh(i4/2a) + D/2a = R^/a^ (12) 

cosh[D/2a)\Jcosh^{D/2a) - 1 + D/2a = R?/a^. (13) 


Defining u = R/a = cosh(Zl/2a) from ([6]) and plugging it in CSI gives u^/vJ — 1 + cosh ^ u = vJ which has a 
solution u ~ 1.211. The value of D/R for which Acatenary = A 2 disks is therefore D/R = 2cosh“^it/u k. 1.055. For 
D/R S (1.055,1.325), the catenary solution has larger area than the two Wilson loops, so it exists only as a local 
minimum. For D/R < 1.055 we evaluate the area (1101) numerically with a from the boundary condition (p|). 
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B. Quadratic Wilson loops 

For quadratic Wilson loops the minimal surface solution would probably look like Fig. [^, j.e., an interpolation of 
Fig. [2]3 and the catenary solution in Fig. [TJa,. We are going to approximate the problem however with the minimal 
area of Fig. and therefore have to minimize the four lateral surfaces. 




FIG. 2. Minimal surface solutions for rectangular Wilson loops. 

This can be done by identifying each of the four lateral surfaces with a catenoid, be., we unroll the circular solution 
in Fig. [T^ and glue four of them together to obtain the minimal surface solution for rectangular Wilson loops shown 
in Fig. [2]3. Starting from the solution (IA8I1 derived in appendix [Aj we can set the integration constant 6 = 0 due to 
symmetry as above, to obtain 


2 /= acosh(x/a) and i?/27r = a cosh(£)/2a), (14) 

where the boundary condition is given by the fact that the perimeter of the catenoid circle at ±Zl/2 (after unrolling) 
has to be equal to the side length R of the quadratic Wilson loop. This changes the range of possible catenary 
solutions and yields an upper bound D/R < 0.88 (IB4I1 . The total surface is now four times the catenoid (fTOjl . be. 

A = 47ra^ [sinh(Il/a) + H/a], (15) 

which equals the Wilson loop areas for D/R = 0.59, hence we again evaluate the minimal area below this value 
numerically with a from the boundary condition (1141) and compare it to the measurements in the next section. 


IV. MEASUREMENTS & RESULTS 

We measure meson-meson potentials in the background of random Z{2) center vortex lines or surfaces in U = 3 and 
U = 4 space-time dimensions. In particular, we measure the correlator of two (flat) Wilson loops with circular and 
quadratic shape, respectively. In fact, these Wilson loops, especially the circular shape, may not represent the perfect 
observable for measuring meson-meson potentials, as creation and annihilation processes of the quark-anti-quark pairs 
can not be neglected. But our main interest lies in the comparison with catenary solutions and therefore these shapes 
are the ideal candidates. 


A. Quadratic Wilson loop correlators in the 4D vortex surface model 

We measure the correlator of two (flat) Wilson loops of size Rx R a.t a distance I? on 16^ lattices of random center 
vortex world-surfaces. The Z(2) vortices contribute a center element (—1) = exp(±i7r) when they pierce a Wilson 
loop. The piercing is unambiguous as the vortices live on the dual lattice, as discussed in section IIIBI The lattice 
string tension of the model is aa^ = 0.755 gives the scale between (minimal) area and the observable. To reduce 
the numerical noise contaminatin g th e measurements as far as possible, the exponential noise reduction technique 
introduced by Liischer and Weisz j55l ] was employed as introduced in for the random center vortex world-surface 
model. It is a multilevel scheme that exploits the locality of the theory by averaging over sub-ensembles, which can 
be applied to the random vortex world-surface model as its action (and the Wilson loop) can be decomposed into 
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sub-lattices. In practice, we use the same number of configurations for the individual averaging steps as detailed 
in to achieve the maximum level of accuracy. However, the correlator runs into the double precision regime and 
breaks down at —log{W{R, R)xW{R, R)x+d) ~ 24, i.e., for individual Wilson loops {W{R,R)) « 10“^^. 

We show 3D plots of the measured correlators and the minimal areas determined from the catenary solution (fTSl) 
or Wilson loop sizes {2R^) in Fig. [3] and 2D cuts of the plots vs. D or i? in Fig. [5j The data shows perfect area law 
behavior and we may see agreement with the catenary solution where the coarse lattice data allows the resolution of 
the effect (Figs. [S)d and d at D = 1 and Fig. [SJ; at D = 1 and 2 resp.) and correspondingly in Figs. [SJd and d before 
the signal is lost in machine precision noise for R and D > 4. 



FIG. 3. a) Z{2) meson (quadratic Wilson loop) correlators —log{W{R,R)xW{R,R)x+D) for various Wilson loop sizes R, or 
distance between quark and anti-quark, and distances D between the meson and anti-meson on 16^ lattices, b) Surface plot of 
the minimum of the two Wilson loop areas 2R? and the catenoid HU. 


B. Circular Wilson loop correlators in the 3D vortex line model 

In the continuous model of random center vortex lines we can define circular Wilson loops, hence we measure the 
correlator of two (flat) Wilson loops of radius i? at a distance D in 16^ volumes with a corresponding string tensions 
aa^ = 0.35 [2l|. The Wilson loops again pick a up a factor (—1) when pierced by a vortex line. In this continuous 
model we can actually go to very small distances where the data precision is still under control, but the catenary 
effect is also very small. We show 3D plots of the data and the minimal areas determined from the catenary solution 
(fTOl) or Wilson loop sizes (27ri?^) in Fig. [3] and 2D cuts of the plots vs. D or i? in Fig. Oand Fig. [6] respectively. The 
results show qualitative agreement with the predicted curves, however they lie significantly below them and can not 
reveal the small difference between cylinder and catenoid surface areas. 



FIG. 4. a) Gircular Wilson loop correlators —log(yV{R)xW{R)x+D) for various Wilson loop radii R and distances D between 
the Wilson loops on 16® lattices, b) Surface plot of the minimum of the two Wilson loop areas 2'kR? and the catenoid (1101) . 























a) D b) 




c) ° d) 




e) D f) 




g) D h) R 

FIG. 5. Z{2) meson (quadratic Wilson loop) correlators —log{W{R, R)xW{R, R)x+d) versus distance R between quark and 
anti-quark (Wilson loop size, right column) and distance D between the meson and anti-meson (left column) on 16^ lattices 
and minimal area solutions (catenary). 
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g) 










FIG. 6. Circular Wilson loop correlators —log(yV{R)xW{R)x+D) for various Wilson loop radii R versus distance D between 
the Wilson loops and minimal area solutions. 
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c) 


e) 





g) 



b) 


d) 


f) 





h) 



FIG. 7. Circular Wilson loop correlators —log{W{R)xW{R)x+D) versus Wilson loop radii R for various distances D between 
the Wilson loops and minimal area solutions. 
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V. CONCLUSIONS 

We measure quadratic and circular Wilson loop correlators in Z(2) center vortex models for the infrared sector 
of Yang-Mills theory, i.e., a hypercubic lattice model of random vortex surfaces and a continuous 2-1-1 dimensional 
model of random vortex lines. We further calculate the catenary solutions for the Wilson loop configurations and 
the corresponding catenoid areas. The measurements show minimal area law behavior and may indicate the catenary 
effects in the hypercubic model, which physically correspond to string surface tension leading to string constriction. 


Appendix A: Beltrami Identity and Catenary Solution 


Starting from the Euler-Lagrange differential equation 


dy dx ^ dyx 


we examine the derivative of / with respect to x 


Solving for the df /dy term gives 


7 o Vx “T rx UxX \ o ■ 

ax oy oyx ox 


r\ yX r\ r\ 7 

oy ox oyx dx 


Multiplying (HD) by yx and substituting the left term with the right hand side of (IA3I) gives 


df df 

dx dx dyx dx dx dyx 


-yx 


dx dy: 

For fx = 0 we derive the Beltrami identity with some integration constant a, 

1 (/ Pi ) — ® ^ f Pi — 

dx oyx oyx 


(Al) 


(A2) 


(A3) 


(A4) 


(AS) 


an identity in calculus of variations discovered in 1868 by Eugenio Beltrami. The quantity / = y^/l + y'^ from Eq.([5|) 
in section UlI Al has in fact fx = 0, so we can use the Beltrami identity to obtain 


iV^ + y''^ -y' 


yy 


\/l + 


72 


= a ^ 2/(1 -b 2 /'^) - yy'"^ = y = ay/l+y’^ 


(A 6 ) 




y = - - 


dx 

dy 


^Jy"^ - a? 


(A7) 


x = a [ — = a cosh ^{-) + b 


y = acosh((^—^), 
a 


(AS) 


which is called a catenary, and the surface generated by rotating it around the x axis is called a catenoid, see Fig. [T] 

Appendix B: Catenary vs. Goldschmidt solution 

To find the maximum value of R/D at which the catenary solutions (lASp resp. ([5]) for the circular Wilson loop 
configuration can be obtained, let p = 1/a, i.e., 


(Bl) 


Rp = cosh{Dp/2). 
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At the maximum value oi D = D* (with corresponding R*) it will be true that dD/dp = 0, hence d/dp of (|B1I) is 

R* = smh{D*p/2){D*/2 + p/2 • dD/dp) = D*/2 ■ smh{D*p/2). (B2) 

Dividing (IBII) at D*, i.e., R*p = cosh(Zl*p/2) by (|B2I) yields D*p/2 = coth(Il*p/2), which has a solution D*p/2 k. 

1.199679 and from (IB2I) we derive the maximum possible value of 

1.32549 (B3) 

for catenary solutions of the circular Wilson loop configuration. For D/R > 1.32549 only Goldschmidt solutions exist. 
For the quadratic Wilson loop configuration the corresponding maximum possible value of D/R is given by 

0.88013. (B4) 
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